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0-^. Abstract 

fj ' We derive a Hamiltonian control theory which can be applied to a 4D symplectic 

C^ \ map that models a ring particle accelerator composed of elements with sextupole 

lyj ■ nonlinearity. The controlled system is designed to exhibit a more regular orbital 

O \ behavior than the uncontrolled one. Using the Smaller Alignement Index (SALI) 

c/3 ■ chaos indicator, we are able to show that the controlled system has a dynamical 

^^' aperture up to 1.7 times larger than the original model. 
^ ! 

MHi Keywords: hamiltonian system, chaos, control, accelerator mappings, 

dynamical aperture 

> 

in 



1. Introduction 

Particle accelerators are technological devices which allow studies at both 
"infinitely small scale" , e.g. particles responsible for elementary forces, and "ex- 
tremely large scale" , e.g. the origin of cosmos. In a simplified approach, such 
devices are composed of basic elements sequence: focusing and defocusing mag- 
nets, accelerating electromagnetic fields and trajectory bending elements as they 
are used in the case of ring accelerators. The resulting dynamics is nonlinear, 
and can be described, in the absence of strong damping, by a conservative 
system. This system can be modelled by a symplectic map built from the 
JH \ composition of several elementary maps corresponding to each basic magnetic 

element. 

One of the main problems found in the dynamics of ring accelerators is 
to study the stability around the nominal orbit, i.e. the circular orbit passing 
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through the centre of the ring. Each component of the ring can be seen as 
a nonhnear map, that deforms the trajectory at large amphtude. Moreover, 
such maps possess stochastic layers whose effect is the reduction of the stability 
domains around the nominal circular orbit (the so-called dynamical aperture 
- DA) [l|. Such behaviors imply that (chaotic) nearby orbits can drift away 
after a few ring turns, eventually colliding with accelerator's boundaries, and 
consequently reduce the beam lifetime and the performance of the accelerator. 

The aim of the present paper is to derive a reliable improvement of the 
stability of the beam by increasing the DA in a simplified accelerator model, 
consisting of only one type of element having a sextupole nonlinearity @,y,|3,lal- 

We work in the framework of the Hamiltonian Control Theory presented 
in [71 y , where two methods to control symplectic maps have been described, 
namely using Lie transformations and generating functions. In the present paper 
we use the former method, that allows direct determination of the new controlled 
map; avoiding the possible problems related to coordinate inversion. 

The aim of control theory is to improve selected features of a given system, 
by slightly modifying its Hamiltonian with the addition of a small control term, 
so that the new system and the unperturbed one are conjugated namely, they 
have the same dynamics. This technique is particularly suitable whenever one 
can directly act on the system and modify it, e.g. in the case of a particle 
accelerator where the addition of a control term in the Hamiltonian function 
can be seen as the introduction of a suitable magnet in the accelerator lattice. 



In our study, we use the Smaller Alignement Index (SALI) [6|, |10|, |11|, [IJ] 
method, which is an efficient indicator for characterising orbits as chaotic or 
regular in Hamiltonian flows and symplectic maps. The SALI is computed 
using the time evolution of two deviation vectors along the studied orbit. 

The paper is organized as follows: we introduce the model in Section [5] and 
present a general result for the control of symplectic maps in Section [3] We 
apply the theory to the symplectic map model of a standard ring accelerator 
in Section U) We briefly recall the SALI chaos indicator in Section [SJ while 
Section [5] presents our numerical results on the behavior of the constructed 
model. Finally, in Section [7] we summarize our conclusions. Further technical 
details can be found in the Appendix. 

2. The model 

Let us consider a system consisting of a charged particle and a simplified 
accelerator ring with linear frequencies (tunes) Qx, Qy, with a localized thin 
sextupole magnet (for more details the interested reader is referred to [2[). The 
magnetic field of this element modifies the orbit once the particle passes through 
it. The basic model is : 
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where xi (xa) denote the deflection from the ideal circular orbit in the hor- 
izontal (vertical) direction before the particle enters the element and X2 (x^) 
are the associated momentum. Primed variables denote positions and momenta 
after the particle left the element. The parameters wi and i02 are related to 
the accelerator's tunes [^ q^ and qy by the relations wi = 2TTqx and W2 = ^nqy. 
The first matrix in ([T|) describes the linear motion of a particle, which corre- 
sponds to a simple rotation in the phase space. The nonlinearity induced by 
the thin sextupole magnet is modeled by the 2"^^ order polynomial expression 
in ([1]). The particle dynamics at the n-th turn, can be described by the se- 
quence {x" ,X2 ,x\P ,x^ )n>Oi where the (n -I- l)-th positions and momenta 
are defined as a function of the n-th ones by ^ . 

The map ([1]) decomposes in an integrable part and a quadratic perturbation, 
respectively the system is associated to the following Hamiltonian (see |Appendix A" ) 



H{xi,X2,x^,Xi) ^ -u)i u!2 and V{xi,X2,xy„XA) ^ — 7r+^i^3- 

n ^^^ 

more precisely ^ can be written in terms of Poisson brackets g as 

x' ^T{x) = e^"^e^^^x. (3) 

Here x — {xi,X2,x^,Xi)'^ , with """ denoting the transpose of a matrix, and by 
definition, for any function / defined in the phase space, {H}f — {H, /} = 

(Vif)'^JV/, with J = \ -, r, ) ' being the symplectic constant matrix, 

e^^V-E^/ ^nd {Hrf = {H}-\{H}f). (4) 

Maps of the form ([T]) have already been studied in [6| where it has been 
shown that chaotic orbits reduce the DA to a hypersphere of radius ~ 0.39 in 
the 4-dimensional phase space (see Figs. 5 and 6 of Q). The goal of the present 
paper is to show that the stability region of the nominal circular orbit can be 
increased once the map (H]) is controlled by an appropriately designed map. 

3. Control theory for symplectic maps 

The aim of Hamiltonian theory is to provide mathematic tools to be able 
to modify the dynamics of a symplectic system. With this, it is possible to 
manipulate intrinsic features, e.g. to reduce the chaotic regions in phase space 
or to build invariant tori. 



^Such parameters have been fixed throughout this work to the values qx = 0.61803 and 
Qy = 0.4152, corresponding to a non-resonant condition (see (Sj). 

■^In the hterature one can sometimes find the alternative equivalent notation {H} = Ljf. 



In the following we will be interested in controlling a quasi-integrable sym- 
plectic map in such a way that it will allow us to obtain a new, controlled map 
"closer" to the integrable part of the original map, and thus increase the stabil- 
ity region around the nominal circular orbit. This is of great importance since 
chaos diffusion channels in phase space with unpredictable consequences in con- 
figuration space. The controlled map is expected to have a smaller number of 
escaping orbits and a larger region occupied by invariant curves in a neighbour- 
hood of the origin. For this purpose, we apply the method presented in [8|, with 
the modification that in the present case the integrable part is not expressed 
in action-angle variables. In particular, the integrable part is a rotation, so the 
present theory applies to perturbations of rotations, instead of maps close to 
identity. 

Let us consider an integrable symplectic map defined through its infinitesimal 
generator |2| H 

a;' = e^^>x, (5) 

and consider the quasi-integrable map perturbation of the former 

x' = T{x) = e^"h^^^x, (6) 

where x G R^^ and 1^ is a perturbation, namely V = o{H). The aim of 
Hamiltonian control theory, is to construct a third map, the control map, whose 
generator F is small (it satisfies F — o{V)). The controlled map 

r,,w=eWe{v}e{n, (7) 

will be conjugated to a map T,, closer to e^^^ than T (see (|lip below). We note 
that the use of the exponential of a Poisson bracket, ensures that such maps are 
symplectic by construction. 

To be more precise, let us define the unperturbed map 

A-' - e-{^> , (8) 

and observe that (1 — ^~^) is not invertible, since its kernel contains any smooth 
function of H. Thus we assume the existence of a "pseudo-inverse" operator, Q, 
that should satisfy (see [8| for details) 

5(1- A-^) g = Q. (9) 

At this point we can define the non-resonant and the resonant operators 

TV := (1 - yl-i) 5 and TZ := I - J\f , (10) 

which are projectors, i.e. Af^ = Af and TZ^ = TZ. 

Our main theoretical result can be stated in the following theorem: 



•^In Q a similar theory has been developed for a general symplectic map devoid of an 
infinitesimal generator. 



Theorem 3.1. Under the above hypotheses and defining S ~ QV we have 

e<^>re*He-^^^-e{^>e{^^>:=T,, (11) 

where 

T,,w = eWe{v}em, (12) 

with a control term given by 

^{F} = ^-{v}^{(J^-g)v]^{nv}^{gv} (^3) 

Remark 3.2 (Warped addition). Let us define as in (^] the warped addition, 
{A} © {B}, of two operators by 

ei^>e{^} := e{^>®^^> . (14) 

An explicit formula can be obtained using the Baker- Campbell-Hausdorff for- 
mula l^], where {A} © {B} is a series whose first terms are 

{A} © {B} = {A} + {B} + 1 {{A}{B} {B}{A}) + ..., (15) 

hence the warped addition is a deformation of the usual addition between oper- 
ators. 

Proof. Using this warped addition, we can rewrite the controlled map into the 
form 

r,,,,=eW®{^}®m, (16) 



where the control term ([13)) becomes 

{F} = ~{v} © {{jv -g)v}® {7^F} © {gv} . (i7) 

From (1101) we have: 



M-g = -A-'g. (18) 

One can easily prove (see appendix A of 8]) that 

h{H}) © i^{S}) © {H} = -{A-'S} , (19) 

hence, recalling the definition of S and (fT8| , we can rewrite (fT7|) as 

{F} = ~{V} © {~{H}) © (-{5}) © {H} © {TZV} © {S} . (20) 

By rearranging the terms we can easily get 

{S} © {H} © {V} © {F} © (-{5}) = {H} © {nV} , (21) 

which is nothing but ()lip rewritten using the warped addition. D 



Remark 3.3. Let us observe that the control term is, as required, small com- 
pared to V . In fact from (J17p and by using the approximated formula for the 
warped addition (|15p . we obtain 



{F} = -{v}(B{{M-g)v}(B{nv}(B{gv} (22) 

= -{v} + {Nv}-{gv} + {nv} + {gv} + o{v)^o{v), (23) 

where we use the relation J\f + TZ = 1 . 

Under the assumption of absence of resonances, i.e. TZV = 0, the o{V) term 
in the control map can be explicitly computed to give (see \Appendix ^ 

F^^{v}gv + o(y^). (24) 

4. The control term for the non-resonant map 

In this section we derive the controhed map presented in section [3] in case 
of maps of the form ([T]). To implement the theory we need to diagonalize the 
operator {i?}. For this reason we introduce complex variables 

Ci = ^2 + ixi and C,2 — Xi + ix^ , (25) 

and rewrite H as 

i/(Ci,C2) = -^|Cip-flC2p. (26) 

The Poisson bracket with H now takes the form 

{H} = lu, (Ci d-^^ - Ci 9c. ) + *^2 (C2 % - C25c. ) . (27) 

Hence, for any n — (rii, 71,2) G N^ and rn = {mi, m2) G N^ we obtain 

{H]Cr = ^ (wiTOi - c^mi + t^2TO2 - ^^2^2) C"C" = ^^ ■ {m - «)("("■ (28) 

Here we introduced the vector uj — (wi,W2) and use the compact notation 
C" = C"X2'^ for the complex vector C — (Ci,C2)- The operator {H} is diagonal 
in these variables and thus map ([8]) is straightforwardly obtained as 

Once we have this map, we can compute the operators g, J\f and TZ. For all 
n and to e N^ \ {0}, such that n ^ rn and 

w • (m - n) 7^ 2kTT Vfc G Z (30) 

(which defines the non-resonance condition), we get 

gcr = z — hj^^Afcr , (31) 

^ ^ 1 Q — iuj-{m—n) ^ ^ ' \ / 



with 



MCC = C '" if '^ • ("^ - ") ?^ 2fc7r Vfc e Z 
" otherwise 



In the rest of the paper, we will assume the above non-resonant condition (|30p 
to hold for the considered values of qx and qy. Let us remark that this is not 
a limitation of the actual theory, but just a working assumption. We could 
equivalently have chosen to work in the resonant regime, using a different control 
term suitable for resonant dynamics. 

The explicit computations are as follows : first we need to express V in terms 
of complex variables; next we compute S — QV, and transform back to the 
original variables. Finally we compute the exponential e^^^\ Actually, even if 
QV is a polynomial of degree three in the variables a;i, 2:2, X3,Xi the map e^^^^ is 
given by an infinite series. The terms of this series can be sequentially computed, 
but the degree of complexity (i.e. the number of involved terms) increase very 
fast. To simplify the computations we decided to use the approximated generator 
of the control map (pi)) already truncated at order 2 

F2 = l{V}gV. (33) 

We note that F2 is composed by about twenty terms. A detailed discussion of 
the whole procedure needed to obtain F2, as well as its explicit formula, are 
presented in [Appendix A[ 

We now face another difficulty, namely the computation of the control map 
from the generator F2. This is equivalent to perform the sum 

e^^^> = l + {F2} + i{^2F + ..., (34) 

whose complexity once again grows very fast. We thus introduce a second 
approximation to our construction, by computing only a finite number of terms 
in the above sum. So, we define a truncated control map of order k 

C,{F2) = Y.^> (35) 

and a truncated controlled map of order k 

n{F2) = e{^>e{^>C,(F2) = rC,(i^2) • (36) 

Since the exact control map e^^^^^ given by ([M)) is symplectic, the controlled 
map ([7]) will also be symplectic. On the other hand, we cannot expect the fc-th 
order control map Ck{F2) to be symplectic. We know that a map is symplectic 
if its Jacobian matrix A verifies (in its definition domain) the equality 

A^JA - J == , (37) 

Thus, in order to check the symplecticity defect of Tk{F2) we compute the norm 
Dk of matrix AjjAk — J, where Ak is the Jacobian of Tk{F2) given by (|M)) . 



The results are presented in Fig. [T]for orders fc = 1 up to fc = 5, in the region 
{xijXs) £ [—1,1] X [—1,1], a;2 = a;4 = 0. The results indicate that Tk{F2) is a 
good approximation of a symplectic map for fc > 4, because we get Dk < 10~* 
for a large portion (> 53 %) of variables values. We note that in the central 
region of the truncated controlled map, where the actual physical process of 
beam's evolution occurs, the symplectic character of the map is established 
even better since there D^ ^ 10~*. As expected, the larger the order fc, the 
closer to symplecticity the approximation is. 
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Figure 1; The simplicity defect of the controlled map Tk{F2) l|36ll . Plot of logj^gDfc, where 
Df^ = ||A^ JAfc — J|| and A^ is the Jacobian of the fc-order controlled map Tk{F2) given by I I36I I. 
for 16000 uniformly distributed values in the square {x\, x^) € [—1, 1] X [—1, 1], 2:2 = 2:4 = 0, 
for (a) fc = 1, (b) fc = 2, (c) A: = 3, (d) k = A and (e) fc = 5. The percentage of orbits with 
logio^fc < -4 is 0.5%, 11%, 28%, 53% and 73% for fc = 1,2,3,4 and 5 respectively. The 



gray scale corresponds to the value of 1 

logj^Q D(j is, and hence the closer the map is to a symplectic one. 



, Di~: the darker the color, the smaller the value of 



The main objective of the addition of a control term is to increase the size 
of the stability region around the central periodic orbit. This increase leads 
to decrease the number of escaping orbitQ as we can see from the results 
presented in Fig. [51 where we plot in black the initial conditions on the square 



"^An orbit {x\ ,x\ ,x), ,x. ] is defined as non-escaping if for all fc < A'' x\ + 

V -^ ^ ^ * J 0<k<N ^ 



^(fc)%^(fc)^ 



(fc)^ 



< R^ for a certain _R (in the simulations we used B? = 10 and N 



10 , 10 ) and escaping otherwise 



{xijXs) £ [—1, 1] X [—1, 1], X2 = X4 = 0, giving rise to orbits that do not escape 
up to 10^ iterations of the map. In particular, we consider in Fig. ^a.) the 
original uncontrolled map ([IJ, and in Figs.[5Jb) to^d), the k order controlled 
map Tk{F2) for k — 1 to k = b respectively. One can easily see that the 
region of non-escaping orbits for the original map is smaller than the one of 
the controlled maps. This observation can be quantified by considering initial 
conditions inside a circle centered at the origin of each panel of Fig. [2] (which 
represent the actual physical plane since the initial momenta are X2 = x\ = Q) 

with radius r = x\ + Xg , and evaluate the number of escaping and non- 
escaping orbits as a function of the circle radius for Tk{F2) with fc = 1 up 
to k = 5. Results reported in Fig. [3] support the previous claim, by clearly 
showing that controlled maps of orders 3, 4 and 5 behave very similarly and 
lead to an increase of the non-escaping region. Let us note that the behavior 
of the controlled maps of orders fc = 1 (Fig. [2][b)) and fc = 2 (Fig. [2jc)) is 
somewhat misleading if it is not analyzed together with the information from 
the symplecticity defect (see Figs. [TJa) and (b) respectively). In fact, these 
maps are strongly dissipative and produce a strong shift of orbits towards the 
origin, preventing them from escaping. This dissipation effect is not physical, 
as it is not observed in real accelerators, and therefore we do not discuss further 
the fc = 1 and fc = 2 controlled maps. 

From the results of Figs. [5] and [3] we see that the addition of even the lower 
order (fc = 3) control term, having an acceptable symplecticity defect, increases 
drastically the size of the region of non-escaping orbits around the central pe- 
riodic orbit. A further increase of the order of the control term results to less 
significant increment of this region, while the computational effort for construct- 
ing the controlled map increases considerably. In fact, Ti(F2) contains around 
100 elementary terms, i.e. monomials in cci . . . 2:4, while this number is almost 
doubled for each order, so that T5{F2) contains around 2000 terms. Also the 
CPU time needed to evolve the orbits increases with the order. For example, 
while the integration of one orbit using Ti{F2) takes about 1.4 times the CPU 
time needed to integrate the original map (H]), the use of T5(F2) needs almost 
21.5 times more. 

Thus we conclude that the T4{F2) controlled map, which can be considered 
quite accurately to be symplectic, is sufficient to get significant increment of the 
percentage of non-escaping orbits, without paying an extreme computational 
cost. 

5. The SALI method 

The Smaller Alignment Index (SALI) [10| has been proved to be an effi- 
ciently simple method to determine the regular or chaotic nature of orbits in 
conservative dynamical systems. Thanks to its properties, it has already been 
successfully distinguished between regular and chaotic motion both, in symplec- 
tic maps and Hamiltonian flows [111 . Il2l llsL [iJ, ll5J . 
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Figure 2: Non-escaping regions of controlled map Tfc (-F2) as a function of the truncation order 
k. 16000 uniformly distributed initial conditions in the square (x\,x^) G [—1,1] X [—1,1], 



(0) 



(0) 
k = Ito k 



- are iterated up to n = 10^ using (a) the uncontrolled map JTJ and (b)-(f) the 
5 order controlled map Tj. (_F2 ) I I36I I , respectively. Initial conditions corresponding 

to non-escaping orbits up to n = 10^ are coloured in black, while escaping orbits are colored 

in white. 
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Figure 3: Percentages of non- escaping orbits for the controlled map Tj.(F2) Il36|l as a function 
of the distance from the origin in the physical space (a'ljXa). We iterate initial conditions in a 



circle of radius r centred at the origin of plane (xi , X3) , with x 



(0) 



.(0) 



: 0, and compute the 



percentages of non-escaping orbits during n 
with k = 1,2, 3, 4, 5, as a function of r. 



10 iterations, for the controlled map Ti^{F2) 



For the sake of completeness, let us briefly recall the definition of the SALI 
and its behavior for regular and chaotic orbits, restricting our attention to 2N- 
dimensional symplectic maps. The interested reader can consult 10] for a more 
detailed description. To compute the SALI of a given orbit of such maps, one 
has to follow the time evolution of the orbit itself and also of two linearly 
independent unitary deviation vectors vl ^v^. ■ The evolution of an orbit of a 
map T is described by the discrete-time equations of the map 



^(«+i) ^y(^(")) 



(38) 



where af^"^ = (a;^" , Xj , ■•■, 2^2'^)'^' represents the orbit's coordinates at the n-th 



iteration. The deviation vectors wj" , v^ at time n are given by the tangent 
map 



^(n+l) AlAn)\ ^ 

v] — A[x^ ) ■ V. 



(n) 



1,2. 



(39) 



where A denotes the Jacobian matrix of map (|38|) , evaluated at the points of 
the orbit under study. Then, according to [lO| the SALI for the given orbit is 
defined as 

(40) 



SALI(n) = min i. 



vl + V2 






}. 



where || • || denotes the usual Euclidean norm and Vi — j^, i = 1,2 arc unitary 
normalised vectors. 
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In the case of chaotic orbits, the deviation vectors vi, V2 eventually become 
aligned in the direction defined by the maximal Lyapunov characteristic expo- 
nent (LCE), and SALI(n) falls exponentially to zero. An analytical study of 
SALFs behavior for chaotic orbits was carried out in [l2| where it was shown 
that 

SALI(n) oc e-'-"'-"^^'', (41) 

with (Ti, cr2 being the two largest LCEs. 

On the other hand, in the case of regular motion the orbit lies on a torus 
and the vectors vi, V2 eventually fall on its tangent space, following a n~^ time 
evolution, having in general different directions. This behavior is due to the fact 
that for regular orbits, the norm of a deviation vector increases linearly in time. 
Thus, the normalization procedure brings about a decrease of the magnitude of 
the coordinates perpendicular to the torus, at a rate proportional to n~^, and 
so vi, V2 eventually fall on the tangent space of the torus. In this case, the SALI 



oscillates about non-zero values (for more details see [ll|). 

The simplicity of SALI's definition, its completely different behavior for reg- 
ular and chaotic orbits and its rapid convergence to zero in the case of chaotic 
motion are the main advantages that make SALI an ideal chaos detection tool. 
Recently a generalization of the SALI, the so-called Generalized Alignment In- 
dex (GALI) has been introduced 16|, [iV | , which uses information of more than 



two deviation vectors from the reference orbit. Since the advantages of GALI 
over SALI become relevant in the case of multi-dimensional systems, in the 
present paper we apply the SALI method for the dynamical study of the 4D 
map H]). 

6. Dynamics of the controlled map 

As already mentioned, the goal of constructing the controlled map Tdri = 
T e i^J is to increase the percentage of regular orbits up to a given (large) number 
of iterations, or equivalently increase the size of the stability region around the 
nominal circular trajectory (i.e. the DA). Because the presence of chaotic regions 
can induce a large drift in the phase space, that eventually could lead to the 
escape of orbits, the achievement of a larger DA can be qualitatively inspected 
by checking via the SALI method the regular or chaotic nature of orbits in a 
neighborhood of the origin (see Fig. 2]). We note that we define an orbit to be 
chaotic whenever SALI(t) < 10~^, and regular for the contrary. 

In Fig. [4] (which should be compared with Fig. 5 of [6]) we observe a strong 
enlargement of the region of regular orbits. This region is characterized by 
large SALI values. In particular, for 10^ iterations of the T4(i^2) map, 54% of 
the considered orbits are regular, while for the uncontrolled map this percentage 
reduces to 33%. This improvement can be also confirmed by visual inspection 
of Fig. [21 where the regions of non-escaping orbits are shown for different orders 
of the controlled map (|36| . 

In Figs.HI^a) and (b) we see that there exist orbits of the T3(i^2) map, which 
are characterized as regular up to n = 10^ iterations, while they show their 
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Figure 4: Stability analysis in the {xi, x:j) plane. 16000 uniformly distributed initial conditions 



in the square {xi,X3) G [—1,1] X [—1, 1], x 



(0) 



.(0) 



: are integrated using the T3{F2) ((a) 



and (b)), and the T4{F2) controlled map ((c) and (d)), up to n = 10** ((a) and (c)) and n = 10^ 
iterations ((b) and (d)). The gray scale represents the value of logj^QSALI for each orbit at 
the end of the integration time. The lighter the color the more stable is the orbit, while white 
color denotes that an orbit escaped before the total number of iterations was reached. The 
black circle indicate the initial condition of the orbit studied in Fig. \5\ 
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chaotic character once they are iterated up to n = 10^. Such orbits correspond 
to the dark regions marked by a black circle in Fig. HJ^b) (for comparison this 
circle is also plotted in all panels of Fig. S]). This discrepancy is absent for 
the Ti{F2) map, which shows almost the same geometrical shape for the non- 
escaping region when we pass from 10'* to 10^ iterations. In order to better 
understand this behavior we followed the evolution of a single orbit with initial 
condition x^^^ = (—0.50, 0, —0.65, 0)^ -inside the black circle in Fig. |4]- for both 
the T^{F2) and the T4{F2) controlled maps, computing the corresponding SALI 
values up to 2 X 10^ iterations. The results are reported in Fig. [5] and clearly 
show that the orbit behaves regularly up to n sa 10^ iterations of the T3(F2) 
map, since its SALI values are different from zero, but later on a sudden decrease 
of SALI to zero denotes the chaotic character of the orbit. This behavior clearly 
implies this is a slightly chaotic, sticky orbit, which remains close to a torus for 
long time intervals {n « 10^), while later on it enters a chaotic region of the 
phase space. It is interesting to note that iterating the same initial condition 
by the Ti{F2) map we get a regular behavior at least up to n = 2 x 10^. 




Figure 5: (Color online) Dynamics of two orbits with the same initial conditions for the 
3'"'* and 4"" order controlled maps. Time evolution of the SALI for the orbit with initial 
conditions x'"' = (—0.50,0,-0.65,0)"^ (see Fig. |4]l, using the T-i,{F2) [(r) red curve] and the 
T4(i<2) controlled map [(b) blue curve]. 

In order to provide additional numerical evidence of the effectiveness of the 
controlled map (j36l) in increasing the DA, we consider initial conditions inside 
a 4D sphere centered at the origin xi = X2 = xa = a;4 = of the map, with 



radius, r^ 



.(0)^ 



.(0)^ 



s^r 



.(") 



We compute the number of regular. 



escaping and chaotic orbits as a function of the sphere radius. The corresponding 
results are reported in Fig.ISJb), while in Fig.ISJa) we reproduce Fig. 6 of ;6l] for 
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comparison. From this figure we observe a strong increase of the DA, since the 
largest sphere containing 100% regular orbits has a radius r « 0.66, while this 
radius was r « 0.39 for the original uncontrolled map. We also observe that 
increasing the total number of iterations from 10^ to 10^ (dashed and solid lines 
in Fig. ini respectively) increases the percentage of chaotic orbits, but the radius 
of the 4D sphere containing only regular orbits does not change significantly. 
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Figure 6: (Color online) Dynamical aperture of the (a) original map {7p and (b) the T/i{F2) 
controlled map \3bt . The percentages of regular [(b) blue curves], escaping [(g) green curves] 
and chaotic [(r) red curves] orbits after n = IC* (dashed curves) and n = 10^ iterations (solid 
curves) for initial conditions in a 4D sphere centred at the origin xi = X2 = x^ = x^ = 0, 
as a function of the sphere radius r. Each point corresponds the average value over 5000 
initial conditions. The largest radius at which the percentage of regular orbits is still 100%, 
is marked by an arrow in each panel. 



7. Conclusions 

In this paper we considered a simple model of a ring particle accelerator 
with sextupole nonlinearity that can be described by a symplectic map. In the 
framework of Hamiltonian control theory, we were able to control the dynamics 
of the original system, by providing a suitable control map, resulting in a small 
"perturbation" of the initial map. This control map has been constructed with 
the aim of DA enlargement of the particle accelerator, and thus improving the 
beam's life-time and the accelerator's performance. 

In particular, the theoretical framework we developed allows a 1-parameter 
family of approximated controlled maps. We performed several numerical simu- 
lations in order to choose "the best" approximated controlled map TkiF^) ([5S)) . 
taking into account the complexity of the map, i.e. the number of terms by 
which it is composed, the CPU time needed to perform the numerical iteration 
of orbits, and the accuracy of the results in terms of the symplectic character of 
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the map. We find that the 4*^* order controlled map T4{F2) is an optimal choice 
for the controlled system. 

Using this controlled map we succeeded in achieving our initially set goal, 
since the T4(i^2) map exhibits a DA with a radius more than 1.7 times larger 
than the one for the original map (see Fig. [5]). 

Appendix A. Computation of the control term. 

The aim of this section is to introduce further details for the construction 
of the control term and of the controlled map, and to provide explicit formulas 
for the interested reader. 

Appendix A.l. Notations 

Let us first introduce some notations and recall some useful relations. 

• Lie brackets and operators. Let X be the vector space of C°° real 
or complex functions of 2N variables {p,q)- For any F, G G X, the Lie 
bracket is given by 

N 

{F, G} := J2 i^p.Fdq^G - d.^Fdp^G] , (A.l) 

where dxif = g^ denotes the partial derivative with respect to the vari- 
able Xi. 

Using the above definition, we can define a linear operator, induced by an 
element i^ of A", acting on X 

{F}: X ^ X 

G ^ {F}G:={F,G} ^ > 

This operator is linear, antisymmetric and verifies the Jacobi identity 

VF,GeX {{F}G} ^ {F}{G} - {G}{F} (A.3) 

• Exponential. We define the exponential of such an operator {F}, by 



oo 



k 



fe=0 

which is also an operator acting on X. The power of an operator is the 
composition : {F}''G = {F}''-^ {{F}G). 

We observe that in the case of the Hamiltonian function H, the exponential 
provides the flow, namely e^^^XQ = x{t), of the Hamilton equations 

jp = -dqH 



q = dpH . 
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• Vector Field. The action of the above defined operators, can be extended 
to vector fields "component by component" 



yF,G,HeX {F} 



G 
H 



{F}G 
{F)H 



(A.6) 



Appendix A. 2. Mappings as time-1 flows. 

We show now that map (H]) can be seen as the time-1 flow of a given Hamil- 
tonian system. More precisely we show that 



f x,\ 

X2 

X3 

\ X4 J 



f cos(a;i) — sin(a;i) 

sin(cL'i) cos(a;i) 





f X, \ 

xi 

x-i 
\xi j 






cos(a;2) 
sin(w2) 



\ 






-sin(w2) 
cos(a;2) j 



( 



xi 



X2 



ei«>e^^> 



X3 
\ X4- 2X1X3 J 



(A.7) 



where 



and 



H{xi,X2,X3,X4) = -OJi- 



-UJ2- 



V{xi,X2,X'i, Xi) = ^ + XiX] 



(A.; 



(A.9) 



Let us observe that H is the sum of two non-interacting harmonic oscillators 
with frequencies oji and a;2, hence its dynamics is explicitly given by 




— Acos{LJit) — Bsa\{ujit) 
~ B cos(wit) + A sin(a;ii) 

— C cos{uj2t) — Dfim{uj2t) 

— D cos{uJ2t) + C sm{uj2t) ■ 



(A.IO) 



By definition y = e^^'x is the solution at time 1 with initial condition 
X— (xi, a;27 2:3, a;4)"^, hence we obtain 



that is 



o{H} 



f x,\ 

X2 

X3 

\X4 J 



= cos(a;i)xi — sin(a;i)a;2 
= sin(wi)a;i -|- cos(a;i)a;2 
~ cos(a;2)a:3 — sin(a;2)a;4 

= Sm{LJ2)x3 + C0S{UJ2)X4 . 



/ cos(wi) — sin(wi) 

sin(wi) cos(a;i) 



\ 



(A.ll) 












X2 


cos(a;2) 


— sin(aj2) 




X3 


sin(a;2) 


C0S(W2) J 




\X4 J 



(A.12) 
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From (jA.gp and the definition (|A.2p we easily get 

{V} := d^^ Vd^, - d^, Vd^^ + d^^ Vd.^^ - d^^ Vd.^^ 
= (xl - xl)dx2 - 2xiX3dxi- 



(A.13) 



This means that once apphed to a vector x only the second and fourth compo- 
nents of {T^}a? will be non-zero and moreover they only depend on the first and 
third components of x, hence {F}^af — 0. We can thus conclude that Vfc > 2 
and Vx € R*, we get {V}''x — 0. Finally using the definition (|A.4p we obtain 



e^^>f = ^ 



{vy 



( 



fc=0 



■X = Ix + {V}x 



\ 



Xi 

X2+ xf- xl 

\ X4- 2xiX3 J 



(A.14) 



Appendix A. 3. Computation of the generator F under the assumption TZV = 0. 

Let us recall that the composition of maps expressed by exponential defines 
the warped addition 

(A.15) 



e{^}e{s} ;:=e^^>®^^>. 



whose first terms are 



{A} ® {B} = {A} + {B} + -{{A}{B} {B}{A}) 



(A.16) 



Using the warped addition with equation (1131) of Theorem 13.11 we obtain 

e{F} = e"'f^>e^(^"^^^^e^'^^^e^^^> = ^-{v}®{{N--g)v}®{TZV}(s{gv} ^ (A.17) 



and thus 
{F} = 



-{V} © {(A/" - g)v} © {TZV} © {gv} = --{v} © {(i - g)v} © {gv} + o{v^) 



,^i{v}{gv}~{gv}{v})~{gv} 



\m 



{{v}gv} 
1 



{gv} + oiv^) 



l{{v}gv} + i {{{v}gv}gv} +o{v-} = -{{v}gv} + o{V') , 



(A.18) 



where we explicitly used the assumption TZV = to remove the third term on 
the the right hand side on the first equation and hence to write AfV = V. We 
are thus able to define the non- resonant control term, up to order V^, to be 



F2^^{v}gv = ^{v,gv} 



(A.19) 
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d 

dxi 

d 


dCi d dCi d 
dxi aci dxi dCi 
dCi d ^dCi d 


. d . d 
d d 

1 


dX2 


dx2 9Ci dx2 dC,i 


5Ci aci 



Appendix A. 4- The operator Q 

To get the explicit formula for F2 we need to compute the expression of Q. 
^From definition ^ the operator Q should satisfy 

g(l-e-t^>)g = a. (A.20) 
To construct it, it will be more convenient to use complex variables 

Ci = 2^2 + ixi and (2=^4 + ix^ . (A. 21) 
Then the function H becomes 

i?(Ci,C2) = -YCiCi-YC2C2, (A.22) 

and using 

a Pi/-. Ft Ft?, a a a 

(A.23) 

(A.24) 

'^2 C(,l 1^1,1 Ui,l 

for (xi,a;2), the operator {H} becomes 
d,,Hd,, - d^.Hd,, = 2i (%//% - %i?9ci) = i^i (Ci% - Ci^cO , (A.25) 

with a similar expression holding for (a;3,a;4). Hence for any n — (ni,ri2) G N^ 
and m = (7711,7712) G N^ we obtain 

{if}C"C = i (wi"^i - ^1^1 + ^2^2 - W2n2) C"C"' = i^-(m-n)C"C , (A.26) 

where we introduced the vector uj = (101,102) and we used the compact notation 
^" = CrC: for the complex vector C = (Ci, C2). 

We note that from the knowledge of the operators' action on such monomials 
^"^™, we can reconstruct the operator action on any regular function. The 
operators are linear and they will be applied on polynomials in the x variable, 
which are nothing more than polynomials in the complex variables. 

Let us now compute the time-1 flow of {H} by using complex variables. 
Starting from (JA.26P and then proceeding by induction, we can easily prove 
that for all fc e N 

{H}''CC = (iw • (m - n)f C"C , (A.27) 

and finally 



{H}^n^m _ \ ^ U^ ' V^^ ^77 ^n^m _ ia-{rh-n) /-n /-m I \ r^Q\ 



k=0 

Similarly e-'f^>C"C™ = e~''^-(™-")C"C'^- 
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Assuming a non-resonance condition 

uj-{7n-n)^ 2kTT Vn 7^ m £ N^ \ {0} and Vfc e Z , (A.29) 

a possible choice for the operator Q is the fohowing 

SCr := ^_^_i(^-.) -A/-Cr (A.30) 



with 



otherwise 



It is easy to check that the operator Q defined by (jA.30|) verifies (jA.201) . In 
order to do so we introduce the compact notation 

Zn,m:=CC and n„^„, := e-'^'^™"") . (A.32) 

Developing the left hand side of (|A.20p and using the linearity of all operators, 
we get 

g(i-e-^"Agz^^„. = (g-Ge-^"^)- — ^ — z„,„ (A.aa) 

-gZn,rn - ^ ^e-^^J Z„^,„(A.34) 



1 — n ' 1 — n 

Zn,m ^ -. ■p. y (nn,m-^n,W-35) 



D 

(i~a 



2 '^n,m /-. |— I \2 ^n^m (^A.duj 



!-□. 



gz^^^ . (A.37) 



Appendix A . 5. Expression of the control term F2 

The function F2 is defined by (jA.19p . where F is a known function. The 
term gV will be computed starting from the previously obtained expression of 
g. To construct F2 we first have to express V in the complex variables (jA.2ip : 

^(Ci, C2) = -^*C? + Y^iCl - Y^iQi + Y^Ql + ^<iCl - ^<iC2C2 

+\KiCI - ^*CiCl + ^*CiC2C2 - ^*CiC2 • (A.38) 

By the linearity of the operator, and by using (|A.30p . we easily compute gV . 
In particular we apply g to each term of (jA.381) . Then using the inverse change 
of coordinates 

xi = xi(Ci - Ci) and X2 = -(Ci + Ci) (A.39) 
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(similar expressions hold for (2:3, a;4) and (C2, Cs)), we can go back to the original 
variables x. The obtained expression after some algebraic simplifications is 

gV = -1/6 csc(3/2wi)[a;2 cos(l/2u;i) + a;i sin(l/2wi)] [a;!^ -3a;32 + a;2^ 
+ (2a:;i —62:3 +2:2 ) cos (cji) — xi a:;2 sin (cji)] 
sin(cj2) 



-1/4 



cos ((x>i — W2) — COS (CJ2) 



_ , (-X2 X4^ + 2:2 2:3^ + 2 xi X4 a::3) sin (1^2) , . . 

cos (CJ2) — cos (wi + (x)2) 

Then the explicit expression of the control term F2 is 

F2 = l/2(a;?-a;^)| -l/6csc(3/2a;i)cos(l/2cji)(a;? -Sx^ + a;^ 

+ (2.Tj^ — 6.T3 + X2) cos(u;i) — a:;ia;2 sin(u;i)) 

— 1/6 csc(3/2a;i) I a:;2 cos(l/2a;i) + xi sin(l/2a;i) J (2x2 + 2x2 cos(a;i) — xi sin((x)i)) 

^u (a^i -a:4)sin(a;2) \ ^ ^^^ xia:3(2a:iX3 - 2x2X4) sin(a;2) (k4,r) 

cos(w2) — cos(a;i + W2) / cos(a;2) — cos(a;i + L02) 
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